Finite-Size Scaling and Universality of the Thermal Resistivity of Liquid ^He near T\ 
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We present measurements of the thermal resistivity p{t, P, L) near the superfluid transition of *He 
at saturated vapor pressure and confined in cylindrical geometries with radii L = 0.5 and 1.0 pm 
(t = T/T\{P) ~ 1). For L — l.O^m measurements at six pressures P are presented. At and above 



Ta the data are consistent with a universal scaling function F = (Lj^SiT^'' [R/Ro), X — {L/f^oY'" t 
valid for all P {Rq and x are the pressure-dependent amplitude and effective exponent of the bulk 
resistivity R and ^ = ^o*"" is the correlation length). 



The modern theory of critical phenomena 0] predicts 
that continuous phase transitions belong to distinct uni- 
versality classes which are determined by such general 
properties of the system as the number of degrees of free- 
dom of the order parameter and the spatial dimensional- 
ity. Within a given class, exponents and amplitude ratios 
are identical (i.e. universal) for all members and indepen- 
dent of irrelevant variables. An example of an irrelevant 
variable is the pressure P of a liquid helium sample at 
which measurements near the superfluid transition tem- 
perature T\ are made. Within a given universality class, 
the dependence of many properties upon certain param- 
eters can be represented by scaling functions which are 
the same for all systems. The present paper is an ex- 
perimental study of the scaling function which describes 
the effect of confinement in a cylindrical geometry with 
radius L on critical properties. For static properties this 
finite-size scaling has been studied by a number of pre- 
cise experiments. For example, the heat capacity near 
the superfluid transition of "^He has been measured for 
confinement sizes which vary by a factor of 500, and the 
data to a large extent can be collapsed upon a unique 
function when properly reduced |^ ^. Even for static 
properties, however, we are not aware of any measure- 
ments which test the universality of finite-size scaling. 
For transport properties there are, to our knowledge, no 
prior experiments which test scaling and universality for 
finite-size effects. There has been only one experiment on 
the effects of confinement on a transport property ^, 
namely the measurements of the thermal conductivity A 
of helium near Tx in cylindrical tubes. These measure- 
ments can be used to derive a scaling function which 
would be expected to be universal, but since they were 
performed only for the one value L = 1/im and only at 
saturated vapor pressure (SVP), they provided a test of 
neither finite-size scaling nor of the universality of the 
derived function. 

We present experimental results for the thermal resis- 
tivity p(t,P,L) = l/\{t,P,L) near Tx{P) of liquid ^He 
confined in cylinders of two different radii and at vari- 
ous pressures as a function of the reduced temperature 



t = T/T\ — l. The use of two confinement sizes allows us 
to directly test finite-size scaling, while the use of differ- 
ent pressures for one size provides a test of universality. 
For bulk helium p{t, P, oo) depends strongly on pressure 
so that a comparison of an appropriate scaling func- 
tion for pit, P, L) at different pressures provides a sensi- 
tive test of universality. These two aspects were tested 
in separate experiments: measurements as a function of 
L were taken at SVP, and measurements as a function of 
P were taken at a single confinement size L = 1.0/im. 

Theoretical predictions for A are still quite limited. 
Monte Carlo calculations give the shape of a scaling 
function, but only to within a multiplicative factor [Q. 
Within its precision this shape agrees well with the mea- 
surements of Ref. 0. Very recently, a one- loop renor- 
malization group (RG) calculation of \{t, P, L) for i > 
and at SVP was carried out by Topler and Dohm , but 
at present there are no such calculations for t < and 
for higher pressures. Thus, in order to provide a broader 
framework for the analysis of our data, we use a phe- 
nomenological approach. We assume that the tempera- 
ture and size dependence of p are separable and that the 
size dependence is a function only of P/^ where ^ = " 
is the correlation length: p{t,P,L) — p{t,P,oo)F [L/^). 
Since p{t, P, oo) goes to zero as t does while p{t, P, L) re- 
mains finite, F diverges at t = 0. To avoid this difficulty, 
we re-define the scaling function as F{X) = [L/S)^^'^ F 
which avoids the divergence at i = 0. Consistent with ex- 
periment 1^ , we have written p for bulk helium as a power 
law p{t, P, oo) — p^t^ with effective exponents x{P) and 
amplitudes po{P). We now have 
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Note that to is the temperature at which the correlation 
length grows to the size of the container, i.e. ^ ~ L. The 
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correlation length has a pressure-dependent amplitude 
^o(-P) and a universal exponent v. The values of £,q,Po, 
and X are known from bulk measurements |^ and are 
summarized in Table 1. 

We used two different thermal conductivity cells. One 
(Cell I) was described in detail elsewhere [Q . It was used 
for measurements of the resistivity as function of P at 
L — 1.0/im. It consisted of two cylindrical metal plates 
made of OFHC (Oxygen-Free High-Conductivity) cop- 
per separated by a stainless steel sidewalk A glass mi- 
crochannel plate (MCP) was epoxied to the inside of the 
sidewall, so that when assembled the liquid helium be- 
tween the plates would be confined to the channels with 
little extraneous liquid between the endplates and the 
glass. A small bulk thermal conductivity cell, called a 
"lambda device" , was attached to the bottom (hot) cop- 
per plate for the determination of Tx in bulk helium. The 
bottom of the lambda device was 1.25 cm below that of 
the confinement cell, and as a result the value of Tx had 
to be corrected for the hydrostatic pressure difference be- 
tween the bottom of the lambda device and the middle of 
the MCP §, p. Since the lambda device was attached 
to the bottom plate, it could only be used before and af- 
ter a data acquisition sequence because the heat applied 
to it necessarily flowed through the confining cell. The 
cell was filled through an overflow volume located on the 
top (cold) copper plate. 

Cell II, used for measurements of the resistivity as 
a function of L at SVP, was designed for use with mi- 
crochannel plates which were surrounded by a solid glass 
ring. Whereas the MCP in Cell I was epoxied into a 
stainless steel sidewall which in turn was sealed to the 
copper endplates with indium gaskets, the glass ring in 
the second type was directly sealed to the copper using 
indium. A stainless steel sidewall was used as a spacer, 
but its length was chosen so that the different thermal ex- 
pansion coefficients of the copper endplates and the stain- 
less steel compensated each other. As a result the force 
applied to the microchannel plate was constant. The 
cryogenic apparatus used with this cell design could ac- 
commodate three thermal conductivity cells, all of which 
were suspended from a common temperature-regulated 
platform. One of these cells was a bulk conductivity cell 
constructed with an open glass ring. It served to locate 
Tx of the bulk fluid. The vertical centers of the cells were 
nearly the same, so that the gravity correction mentioned 
for the first cell was greatly reduced. The three cells were 
independent of each other, so that the bulk thermal con- 
ductivity could be measured at the same time as the ther- 
mal conductivity for two different confinement sizes. The 
fill line entered the bottom of the cell, and the portion of 
the fill line located in the bottom plate was packed with 
0.05 /im alumina powder to suppress the superfluid tran- 
sition. Thus the liquid helium contained in the bottom 
plate was always normal. The fill line was connected to 
an overflow volume on the shield stage, which was main- 



tained a few mK above T\. 

Saturated vapor pressure was maintained in both sets 
of experiments by partially filling the overflow volumes. 
Pressures other than SVP were reached in Cell I us- 
ing a "hot volume" a separate thermal stage filled 
with fluid whose temperature was controlled to regulate 
the pressure in the thermal conductivity cell. The pres- 
sure was measured using a capacitative strain gage 
mounted on the top of the cell. The fluid in the cell and 
the hot volume was isolated from the rest of the cryostat 
by a normally-closed low-temperature valve. 

The resistivity p was computed from the temperature 
difference AT which was measured across the fluid layer 
when a power Q was applied to the bottom plate of 

a cell: p = {A/d) {Q/AT - CwV^ - Rb , where A is 
the cross-sectional area of the fluid, d is the spacing be- 
tween the plates, Cw is the parallel conductance of the 
stainless steel sidewall and the glass of the MCP, and 
Rb is the boundary resistance between the copper end- 
plates and the fluid. The boundary resistance Rb was 
measured far below Tx, where the resistance of the fluid 
layer can be neglected. The size of the correction is rel- 
atively small, and its temperature dependence near Tx 
JlSf was neglected. The parameters A/d and Cw were 
obtained by fitting the measured A(t, P, L) to the known 
A(t, P, oo) several mK above Tx where the effects of con- 
finement are negligible. For Cell I, the pressures at which 
measurements were made were chosen to match those for 
which prior measurements for bulk helium were avail- 
able the value d/A = 0.386 so obtained was found 
to be independent of pressure and agreed with the value 
0.39 previously determined for this cell 0. The 0.5 and 
1 /im data, taken with Cell II, yielded d/A = 0.0781 and 
0.0605 respectively. All values for d/A are in good agree- 
ment with values from gas flow-impedance measurements 
on and electron micrographs of the microchannel plates. 
The values for Cw for each size and pressure are shown 
in Table || (for P — 11.25 bar, there were no bulk con- 
ductivity data, and Cw was obtained by interpolation 
between other pressures). Each conductivity data point 
was assigned to the mean temperature T — T^Qp + Ar/2, 
and a corresponding curvature correction ^ was applied 
to correct for the use of a finite Q and AT. 

The resistivity at SVP is plotted versus t in Fig. |l| for 
two different values of L. The data show the effect of 
confinement, with the smallest size showing the greatest 
rounding of the transition and the greatest increase of 
p{t — 0). The scaling variable F (Eq.Jl^) is plotted versus 
X (Eq. 1^) for the two sizes in Fig. Except perhaps 
for X < —2, the data collapse onto a single curve, thus 
supporting the concept of finite-size scaling. It is difficult 
to tell whether the small difference in F between the two 
data sets below A" ~ —2 is significant or due to unknown 
systematic experimental errors. 

Figure ^ shows p{t, P, L) as a function of t for six dif- 




FIG. 1: Thermal resistivity versus reduced temperature at 
SVP for L — 0.5 fim (open circles) and 1.0/im (open squares). 
The plusses are bulk measurements (Ref. Q) and the dashed 
curve is a powerlaw fit to the bulk data. 



FIG. 3: Thermal resistivity versus reduced temperature 
for L — l.Ofim at SVP (open circles); 6.95bars (open 
squares); 11.25bars (triangles); 14.73bars (inverted triangles) ; 
22.31bars (diamonds); 28.00bars (crosses). The reduced tem- 
perature for each pressure is defined relative to Tx{P). 
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FIG. 2: Scaling function F versus scaling variable X at SVP 
for L = 0.5^m (open circles) and l.O^m (open squares). 



ferent values of P and L = 1.0/xm. The resistivity does 
not drop to zero below t = 0, as is the case for the bulk 
fluid Q. The value of p{t = 0,P,L) varies by nearly a 
factor of three for the pressures used. 

In Fig. ^ the function F is plotted versus X for six 
different pressures. Within our resolution the data col- 
lapse on the same curve, suggesting that a single scaling 
function describes all six pressures. The collapse occurs 
despite the large variation of p at constant t. The values 
of F at X = are given in Table |. 

The thermal conductivity A is plotted on a logarithmic 
scale versus t on a linear scale in Fig. ^ for temperatures 
below Tx{P). It is consistent with an exponential growth 
below Tx{P) as noted previously Q. The amplitudes and 
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FIG. 4: Scaling function F versus scaling variable X 
for L = 1.0/im at SVP (open circles); 6.95bars (open 
squares); 11.25bars (triangles); 14.73bars (inverted triangles) ; 
22.31bars (diamonds); 28.00bars (crosses). 



arguments of the exponential are approximately the same 
for all pressures. These results suggest that A, rather 
than the scaling function F ^ is independent of pressure 
in this temperature region, and that universality breaks 
down below T\(P). 

Aside from testing scaling and universality, an impor- 
tant issue is to what extent detailed theoretical calcu- 
lations can reproduce the conductivity. As discussed 
above, the theoretical information is limited. Monte 
Carlo calculations, which give the shape of the scaling 
function quite well, involve as yet undetermnined param- 
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son, and M. Larson, Phys. Rev. Lett. 84, 4894 (2000). 
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FIG. 5: The thermal conductivity below T\ on a loga- 
rithmic scale versus reduced temperature on a linear scale 
for L = 1.0/im at SVP (open circles); 6.95bars (open 
squares); 11.25bars (triangles); 14.73bars (inverted triangles); 
22.31bars (diamonds); and 28.00bars (crosses). The reduced 
temperature for each pressure is defined relative to Tx{P). 
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FIG. 6: Thermal conductivity X{t = 0) vs L"^ on logarithmic 
scales. The plus is the SVP measurement from Cell I. Open 
squares are data from cell II. The dashed straight line is the 
prediction based on Eq. |l|. The solid curve is the prediction 
by Topler and Dohm 



eters. However, the recent renormalization group calcu- 
lations have yielded results for A at SVP In Fig. ^ 
we show data for X{t = 0) as a function of on log- 
arithmic scales. The phenomenological scaling function 
Eq. predicts X(t = 0) a L-""/" which, for x/iy = 0.656 
is shown by the dashed straight line. The RG prediction 
is given by the solid line(note that the horizontal axis dif- 
fers from the original; here, L is the radius of the channel, 
not the diameter). Although it comes modestly close to 
the data points, it does not follow a pure power law, as 
manifested by curvature of the solid line in the figure. 
Near L ~ 1 the effective exponent (i.e. the local slope 
of the line in Fig. ||) is close to 0.56, which differs sig- 
nificantly from the prediction based on Eq. |l|. Although 
the data points tend to favor the larger exponent x/v, 
they do not cover a sufficient range of L to be decisive. 
Future plans call for the measurement of the resistivity 
at larger L, with the largest (50/im) to be flown on the 
International Space Station, extending the range covered 
to two decades in L. 
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